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SUMMARY 

The rotordynamlc behavior of turbomachinery Is critically dependent on 
fluid dynamic rotor forces developed by various types of seals and bearings. 

The occurrence of self-excited vibrations often depends on the rotor speed and 
load. Misalignment and rotor wobbling motion associated with differential 
clearance have often been attributed to stability problems. In general, the 
rotative character of the flowfield is a complex three-dimensional system with 
secondary flow patterns that significantly alter the average fluid circumferen- 
tial velocity. A multidimensional, nonorthogonal , body-fitted-grid fluid flow 
model is presented that describes the fluid dynamic forces and the secondary 
flow pattern development in seals and bearings. Several numerical experiments 
were carried out to demonstrate the characteristics of this complex flowfield. 
Analyses were performed by solving a conservation form of the three-dimensional 
Navier-Stokes equations transformed to those for a rotating observer and using 
the general-purpose computer code PHOENICS with the assumptions that the rotor 
orbit is circular (nonzero dynamic eccentricity) and that static eccentricity 
is zero. These assumptions have enabled a precise steady-state analysis to be 
used. 


Fluid injection from ports near the seal or bearing center increased 
fluid-film direct dynamic stiffness (Kp = K - Mu 2 (l - X ± X e )) and, in some 
cases, significantly increased quadrature dynamic stiffness (Kq = Dco(l - X ± 
\ q ) - K X y). Under certain conditions Kg was less than zero, implying insta- 
bility (necessary but not sufficient condition) and geometric configuration 
dependency (refs. 1 and 2). Injection angle and velocity could be used for 
active rotordynamic control; for example, injection, when compared with no 
injection, increased direct dynamic stiffness, which is an important factor 
for hydrostatic bearings. 



Three turbulence models were tested: Prandtl mixing, standard k-e, and 

modified k-e. The Prandtl model underpredicted and the standard k-e model 
overpredicted secondary flow zone. The modified k-e model was postulated to 
provide better predictions of secondary flow zones and dynamic stiffness 
(direct and quadrature). 

The Richardson parameter related turbulence production and dissipation 
within the passage; positive Richardson numbers were associated with zones of 
secondary flow and represented regions where turbulence dissipation exceeded 
production, an important observation for stable seal and bearing designs. 
Designs should provide low production near the stator. 


INTRODUCTION 

The characterization of turbulence in high shear flows is at best diffi- 
cult, perhaps impossible. Such flows can and do occur within the passages of 
rotating cylindrical configurations with or without axial pressure drop. For 
industrial applications these passages represent seals and bearings, respec- 
tively, but could equally well represent high-speed propulsion systems or bal- 
listic projectiles. A major difference between characterizing high-shear flows 
in rotating systems and those of an airframe, for example, is the development 
of centrifugal and Coriolis forces and the coupling between an unsteady energy 
source and the working fluid. 

In turbulence modeling of high-shear flows (seals or bearings) researchers 
assume the validity of the universal friction laws, which are based on low- 
shear-flow data. Such modeling ranges from indirect inference of irreversible 
momentum loss to production and dissipation mechanisms. For example, Bentley 
et al . (ref. 3) use a rotordynamic model with a bulk flow parameter to charac- 
terize circumferential velocity (indirectly related to wall or fluid shear); 
Nelson and Nguyen (ref. 4) use gas dynamics equations with Hirs' bulk model 
for friction; Tam et al . (ref. 5) use Prandtl's model; Nordman et al . (ref. 6) 
use the standard k-e model. Each researcher obtained good agreement with the 
experimental data of Childs (ref. 7). Such agreement is necessary but not suf- 
ficient to justify a model or to characterize the turbulence of high shear 
flows . 

This paper compares three models, Prandtl, standard k-e, and modified 
k-e, and relates the results, along with fluid injection, to seal and bearing 
dynamics. The characterization of turbulence is not established as the Navier- 
Stokes solver uses wall functions, derived from the universal law, and coarse 
grid in the near-wall region (NX x NY x NZ = 12 x 6 x 16). 


SYMBOLS 

C(j empirical constant 

D damping coefficient, N s/m 

E constant (integrated function of surface roughness) 

K[) fluid-film direct dynamic stiffness, K - Mw 2 (l - X ± \ e ), N/m 
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fluid-film quadrature dynamic stiffness, Dud — X ± Xg ) — Kyy , N/m 

cross-coupled stiffness, N/m 

turbulence production parameter, (m/s) 2 

mass, kg 

node spacing 

production of turbulence kinetic energy, eq. (4), kg/m s 3 
Richardson parameter 

turbulence dissipation source term, kg/m s 4 
production source terms in two-equation model, kg/m s 3 
dissipation source terms in two-equation model, kg/m s 4 
time, s 
velocity, m/s 

circumferential velocity, m/s 
average velocity, m/s 
coordinates, m 

turbulence dissipation parameter, m 2 /s 3 
von Karman constant 
average velocity parameter 

injection or preswirl average velocity parameter 
viscosity, kg/m s or Pa s 
density, kg/m 3 
production Prandtl number 

dissipation Prandtl number 

wall shear stress, Pa 
perturbation angular speed, rad/s 

rotational angular speed, rad/s 
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TURBULENCE MODELS USED IN SIMPLE TURBULENT FLOWS 


Turbulent flows are generally modeled by using a two-equation k-e turbu- 
lence model (refs. 8 and 9). In the k-e model the local state of turbulence 
is characterized by two scalar quantities, viz, turbulence kinetic energy k 
and its dissipation rate e. The local turbulent eddy viscosity is com- 
puted from k and e as follows: 




t = 



( 1 ) 


where p is the fluid density and Cq an empirical constant (ref. 8 ), all 
derived from flat-plate data assumed to be in local equilibrium. Both k and 
e are here supposed, as in the original publications (refs. 8 and 9), to be 
governed by conventional "transport" equations: 


Dk 3 (V t 9k_^ 

p Dt ' U 3 v 


* P k - pc 


De_ JL Bt_\ 

P Dt " 9x i l a e 8x J 


+ k (C l P k C 2 pe) 


P k = p t 




( 2 ) 


(3) 


(4) 


The empirical constants used in the calculations are those recommended in the 
original references (refs. 8 and 9), Ci = 1.44, C 2 = 1.92, Cq = 0.09, 

< 7 ^ = 1.0, and a e = 1.3. This model has been applied to seals with good 

results (ref. 4). 


Boundary Conditions at Wall 


For flows near walls it can be presumed (1) that the velocity profile 
obeys the universal logarithmic law for smooth walls and includes a rough-to- 
smooth correlation adjustment for wall roughness, ( 2 ) that the turbulence 
energy is proportional to the wall shear stress, and (3) that the length scale 
is proportional to the distance from the wall. 


u 





k = 


x/p 

VTi 


e = 



(5) 

( 6 ) 


(7) 
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where the wall shear stress r is determined from the velocity u in the 
near-wall region and k = 0.435. 


Influence of Curvatures and Swirls on Turbulent Flows 

Experimental investigations of flows with significant streamline curva- 
tures, swirling flows in particular, reveal that the extra centrifugal and 
Coriolis forces can have a surprising effect on turbulence structure and mean 
flow patterns (ref. 10). It has been observed that turbulence is generally 
augmented near the concave and reduced near the convex surface of curved chan- 
nels. Figure 1(a) illustrates the typical swirling velocity profile in a 
curved channel and indicates the relation between the gradient du/dr of the 
swirl velocity and the turbulence intensity (i.e., du/dr < 0 and du/dr > 0). 
Reduced outward diffusion of turbulence kinetic energy from the convex inner 
wall, as well as the net force of both the radial pressure gradient and the 
centrifugal force, prevent the formation of large-scale turbulence. The effect 
of the concave wall is just the opposite. 

In forced swirl, such as that formed by a rotating shaft (fig. 1(b)), the 
centrifugal forces near the inner wall support the growth of vortical struc- 
tures even more than in curved-channel flow. The angular velocity gradients 
are generally negative all across the flowpath, with the largest absolute val- 
ues near the stator. It is in that region where the largest turbulence will be 
generated, the "separated" region (fig. 1(b)). This situation is inherently 
unstable since high turbulence levels weaken the secondary flow patterns. How- 
ever, introducing preswirl or multi injection in the seal or bearing configura- 
tion can often cause a counterswirl ing vortex to form near the housing wall. 
This has been known in practice to have a stabilizing effect on seal, bearing, 
and rotor dynamics (ref. 5). In this case the largest angular velocity gradi- 
ent is shifted inward toward the shaft. Within the countervortex zone the 
Richardson parameter and the gradient of angular velocity change sign, reduc- 
ing the turbulence levels in that region and enhancing the system stability. 


The Modified k-e Turbulence Model 

Existing turbulence models (refs. 8 and 11) can successfully predict only 
simple boundary layer or mildly recirculating flows. For coupled swirling 
flows they fail to capture the essential physics of usually anisotropic, nonho- 
mogeneous, and three-dimensional turbulence structures. Reported attempts at 
predicting rotating turbulent flows indicate that complex, full Reynolds stress 
equations provide no better predictions than those obtained with simple correc- 
tions to the k-e model (ref. 11). The general consensus is that the empiri- 
cal constants in semiempirical models (such as standard k-e) are invalid for 
complex flows. 

For this reason the authors' previous publication (ref. 5), where a simple 
Prandtl mixing-length model was used to represent the turbulence in seal and 
bearing flow passages, may not be correct. Herein an attempt was made to eval- 
uate a modified k-e model and compare its results with those of Prandtl 
mixing-length and standard k-e models. 

Several investigators have attempted to extend both algebraic Reynolds 
stress (ARS) and k-e models to turbulent rotating flows. For a comprehensive 
review consult the work of Lakshminarayana (ref. 11). 
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Most of the relevant approaches are based on Bradshaw's suggestion (ref. 10) 
that the streamline curvature effects can be modeled by modifying the turbu- 
lence length scale. Two most relevant approaches are those of Rodi (ref. 9) 
and Launder et al . (ref. 12). Both incorporate the "rotating" correction into 
constants of the turbulence dissipation rate source term, viz, 

S e = ( C l P k " C 2 pe ) (8) 

where is the generation of turbulence kinetic energy k and C] = 1.44 

and C 2 = 1.92 are the "universal" constants. Rodi's extension introduces a 
correction to the model constant Ci in the generation part of the e-equation 
source term: 


C ] = C^l + C'Ri) (9) 

Launder 1 s extension modifies the destruction part of the e-equation source 
term in a similar manner by modifying the C 2 constant as 

C* = C 2 (l + C^Ri) (10) 

where = 0.9 and C£ = 0.2, respectively, and Ri is the Richardson param- 
eter. The models differ in their specifications of the Richardson parameter, 
which they define as 


(7)- 

Both respond to the angular speed u^/r 2 , the swirl velocity gradient, and the 
turbulence time scale, which differs in preceding definitions. 

Both approaches were also recently compared by Chen (ref. 13) for confined 
swirl jet flow, with the conclusion that Launder 1 s correlation gives the best 
prediction in the recirculating zone. This approach was also used in the 
present analysis. 



Rodi (modifies C, ) (11) 

Launder (modifies C 2 ) (12) 


Present Analysis 

Analysis of the correcting factor indicates that in the regions of nega- 
tive velocity gradients, Ri is also negative and smaller e decay intensifies 
the turbulence level. The negative Richardson parameter starts increasing from 
the shaft surface, and its magnitude becomes maximum in the separation region. 
Then the gradient of the Richardson parameter drops sharply and changes to the 
positive Ri parameter, where the turbulence level is damped and a secondary 
recirculating zone appears. 
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Current Development of Turbulence Models 

The comparative verification studies of Chen (ref. 13) and also of Roback 
and Johnson (ref. 14) indicate that the predictions show correct trends but are 
at best qualitative. It is the authors' opinion that further turbulence model- 
ing for swirling flows should follow a more rigorous path rather than modifying 
"universal" turbulence constants. Two visible improvements of the standard 
k-e model have been recently reported and are worth testing for rotating 
flows; these are the extended model of Chen and Kim (ref. 15) and the multiple- 
scale model of Chen (ref. 16). In the first model two time scales are intro- 
duced in the dissipation rate equation: ( 1 ) the production-range time scale 

tp = pk/P^; and (2) the dissipation-range time scale to = k/e. With these 
scales the source term of the e-equation is 


S 

e 


l P k k ‘ C 2 P 


e 

r + 


C 

_s\ 

pk 


or in terms of the time scales 


where 


S 

e 


C,P, 





(13) 


(14) 


t 


D ~ 


k 

e 


The last term represents the energy transfer rate from large-scale turbulence 
to small-scale turbulence and is controlled by the production-range time scale. 
The recommended constants are C] = 1.15, C 2 =1.9, and C 3 = 0.25 (ref. 15). 
Results obtained with this model show significantly improved accuracy for 
swi rl i ng flows . 

The second model, the multiple-scale k-e model, has recently been suc- 
cessfully implemented by Chen (ref. 16). Following Kolmogorov (who postulated 
that the turbulence spectrum comprises independent production, inertial, and 
dissipation ranges), Hanjalic and Launder (ref. 17) have proposed a multiple- 
scale model in which separate transport equations are solved for the turbulence 
energy and dissipation rate across the spectrum. More recently Chen (ref. 13) 
has improved Launder's original formulation by ensuring an overall kinetic 
energy balance within the spectrum. The model includes two sets of k and e 
equations for large-scale energetic eddies (P = production range) and for dis- 
sipative eddies (I = inertial or transfer range), as shown in figure 2 . The 
source terms of the corresponding transport equations for this model are 


kP = P k pe P ) 

(15) 

J Production 


ki ■ p (, p - c i> ) 

(16) 
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( 17 ) 


S eP " <C P 1 P k C P 2 pe P )7t P ) 

\ Dissipation 

S el = p(C I 1 C P " C I 2 e I )/t P | (18) 

where t p = k p /e p is the production time constant and is defined differently 
than in Launder ' s model . 

The coefficients are 


C p = 1.6, 
*1 


C 


kP " a kl 

= 1 


(19) 

= a el = 

1 .22 


(20) 

= 1.8- 

no 

l - kj/kp 

(21) 

u . c ■ 

l + kj/kp 


1.15, C T =1.8— (22) 

l 2 e P 


By using these model equations the Reynolds stresses are expressed according 
to 


P<u! ,uj> 


- p t 


0.5 



3 pk5 ij 


where 


k = kj + k p 


and the eddy viscosity is given by 


= c p P (k p + 


(23) 


(24) 


(25) 


with C p = 0.09- 

Results obtained with this model show significantly improved accuracy over the 
standard k-e model (eqs. (2) and (3)), have good prospects in rotati ng-f low 
applications, and represent a logical extension of this work. 


FLUID INJECTION AND TURBULENCE MODEL RESULTS 

Control of the integral, average circumferential velocity can lead to con- 
trol of instabilities in rotating machines (Benkert and Wachter (ref. 18), 

Miller (ref. 19)); swirl brake and antiswirl applications (Kirk (ref. 20)); 
swirl and injection (Brown and Hart (ref. 21) and Bently and Muszynska (ref. 3)); 
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injection (Hendricks (ref. 22)); friction and wall surface control (Nelson and 
Nguyen (ref. 4) and Von Pragenau (ref. 23)); preswirl (Childs (ref. 7)). 


Fluid Injection 

The fluid injection geometry is given in figure 3. Tam et al .' (ref. 5) 
found that for a dynamic eccentricity of 0.8 with precessional speed equal to 
rotor speed, tangential fluid injection from four jets (at 0°, 90°, 180°, and 
270°) angled at 3:1 (72°) significantly increases both direct, Kg = K - M[« p - 
(X + X 6 )cor] 2 , and quadrature, Kq = DCup - (X + X e )u>R] - K X y, dynamic stiff- 
nesses but that radial injection tends to be unstable (fig. 4). Tam also found 
that above 450 rad/s injection against the angle of rotation (-3:1, or - 12 °) 
produces larger quadrature dynamic stiffness but smaller direct stiffness than 
injection in the direction of rotation. 

In the present study fluid injection against rotation at -1:1 (-45°) first 
produced an increase in dynamic stiffness (to 350 rad/s) and then decreased it 
to significantly less than that for injection with rotation at 1 : 1 (45°). How- 
ever, for either case the dynamic stiffness was larger than Tam (ref. 5) found 
for injection at ±3:1 (±72°). This implies the existence of an optimum injec- 
tion angle. The quadrature dynamic stiffness for injection against rotation 
at -1:1 (-45°) increased significantly with rotor speed, surpassing the other 
injection results. This implies enhanced configuration stability when fluids 
are injected against rotation. Figure 5 provides some qualitative details of 
the flowfield at L/4, L/2 - E, L/2, L/2 + E, and 3L/4, where L is the con- 
figuration length and E = L/20. The rotor speeds selected were 245 and 532 
rad/s (2344 and 5085 rpm) 1 with fluid injection at an angle of 1 : 1 (45°) both 
against and in the direction of rotation (figs. 5(a) to (d), respectively). 
Preswirl was not considered. Note the zones of secondary flow, and in particu- 
lar note the forward recirculation along the axis of rotation; the intensity 
of these zones increased when injection was against rotation and proportional 
to rotor speed. Note the emergence of discontinuous circumferential secondary 
flow zones at the higher rotor speed. It appears that with 1:1 injection the 
circumferential component was sufficiently weak to permit flow reattachment 
but the radial component was strong enough to "block the flow" and force fluid 
to recirculate upstream. Such secondary flows have not been found experimen- 
tally but have some parallel in the pocket bearing visualization work of Braun 
(private communication). 

Comparing figures 4 and 5 indicates that at reduced rotor speeds injection 
against rotation enhanced the secondary flow zone, reducing the integral aver- 
age circumferential velocity and increasing direct and quadrature dynamic 
stiffness. However, at the higher rotor speeds injection against rotation 
decreased direct but increased quadrature dynamic stiffness. This suggests 
that stiffness and damping could be decoupled in designing seals and bearings. 
Although injection angle and mass flow are important, other parameters such as 
clearance, radius, static eccentricity, rotor speed, and fluid properties 
remain to be investigated. 


Uhe working fluid is bromotrifluoromethane (CBrF 3 ), the shaft diameter is 
7.7 cm, and at these rotational speeds the Reynolds numbers simulate SSME class 
seal s. 
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Turbulence Model i ng 


To investigate how turbulence modeling affects dynamic seal flows with 
inlet swirl, three models were used: Prandtl , standard k-e, and modified k-e. 

The steady three-dimensional Navier-Stokes equations using the PHOENICS solver 
(ref. 24) provided numerical results. The configuration and parameters are 
given in figure 6. The effective viscosity at three locations, 90° from the 
maximum clearance, approaching maximum clearance, and approaching minimum 
clearance is shown in figure 7. In most instances the general character of the 
profile was M-shaped. The peaks were prominant for the Prandtl model and 
underpredict the secondary flow zones as shown in figures 8 at L/8 and L/4 for 
a rotor speed of 5085 rpm, no preswirl, and no injection. The peaks were high- 
est in the maximum clearance zone and nearly the same at 90° but lower when 
approaching the minimum clearance, indicating possible relaminarization 2 of the 
flowfield within the minimum clearance region as postulated by Braun et al . 

(ref. 25). 

The standard k-e model also peaked but at a reduced magnitude and tended 
to overpredict secondary flow zones. It was postulated that the modified k-e 
model, incorporating the Richardson parameter (Ri = (k/e ) 2 |(u@/r 2 )| dru@/dr 
(eq. (12)) and producing results bounded by the other two models, would more 
accurately predict secondary flows and effective viscosity. In the absence of 
secondary flows turbulence production exceeded dissipation from the shaft to 
the wall (i.e., the Richardson parameter was negative and decreased, fig. 9(a); 
in some instances it may achieve a magnitude of 5). With secondary flows 
present, the Richardson parameter was first negative and decreasing but near 
the wall it abruptly increased and became positive (fig. 9(b)), indicating that 
turbulent dissipation exceeded production. This was also the region of second- 
ary flow where direct and quadrature dynamic forces enhanced stability, imply- 
ing that regions of positive Richardson parameter improved rotor stability. 

Regions of positive and negative Richardson parameters, so called Richard- 
son flux, are illustrated in figure 10, along with circumferential velocity at 
L/4 from inlet and L/4 from exit for rotor speeds of 3252 and 5085 rpm and 
preswirl with and against rotation; the modified k-e model was used. The 
positive Richardson parameter zones are closely associated with regions of sec- 
ondary flows, with negative zones being void of secondary flows. Note, how- 
ever, that for both 3252 and 5085 rpm with preswirl against rotation, the Rich- 
ardson flux became negative but the secondary flow zone persisted. The strength 
of these secondary flow zones depends among other things on eccentricity and 
could produce negative Richardson parameters at smaller eccentricities. How- 
ever, as the shaft moves toward the housing, positive Richardson parameters 


2 Relaminarization within a flowfield approaching the minimum clearance, in 
particular near the rotor, would significantly reduce both the viscous forces 
and their perturbations within the minimum clearance, implying enhanced stabil- 
ity: The turbulent flowfield outside this zone, with viscous forces an order 
of.magnitude or more larger, would be rebuffed or turned around, as in second- 
ary flow that numerically has been shown to occur. A superlaminar flow could 
also be engendered in the same manner, and the flowfield would appear to be 
similar. Downstream of the minimum clearance the flow would then be character- 
ized by transition from laminar to turbulent flows with flow separations or 
cavitations with either turbulent or superlaminar flows. 
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could be produced (i.e., first turbulent production exceeds damping, then tur- 
bulent damping exceeds production). The respective forces would tend to be 
first destabilizing and then stabilizing, causing drift followed by an abrupt 
movement of the shaft toward the centered position (i.e., an instability 
quench). This type of oscillation is apparently commonplace in turbomachinery 
and was pointed out to the authors by Childs (ref. 26). 

It is interesting to note that the direct dynamic stiffness forces pre- 
dicted by the standard k-e model were greater than those predicted by the 
modified k-e model for preswirl against rotation (fig. 11(a)). However, 
there was practically no difference in the quadrature dynamic forces 
(fig. 11(b)). At this time, we have no explanation for this result. 

In comparing direct and quadrature dynamic forces, for the modified k-e 
model, the direct dynamic stiffness was less and had a larger negative slope 
for preswirl against rotation than for preswirl with rotation (fig. 12(a)). 

But the quadrature dynamic stiffness was greater and had a larger positive 
slope for preswirl against rotation than for preswirl with rotation (fig. 12(b)). 
With the formalization of the measurements of Morrison et al. (refs. 27 and 
28), incorporation of the two-region model, and comparison to other Navier- 
Stokes solvers, the characterization of turbulence may be possible. 


SUMMARY OF RESULTS 

In this paper some effects of fluid injection into high-shear flows are 
discussed, and three turbulent models used to assess high-shear flows such as 
those in seals or bearings are examined. 

Injection at 1:1 (45°) engendered a recirculating zone upstream of the 
injection plane and increased direct dynamic stiffness (magnitude and slope) 
over that for injection at 3:1 (72°). For injection at 1:1 the quadrature 
dynamic stiffness was negative at low speeds (unstable) and became positive at 
high speeds with a rather large slope for injection against rotation. These 
results suggest that designing for stiffness or damping may be configuration 
specific; further, variable injection angle and velocity could become the 
basis of a practical dynamic control system (i.e., the increase in stiffness 
is important to hydrostatic bearings). 

Of the three models tested (Prandtl, standard k-e, and modified k-e) the 
Prandtl model tended to underpredict and the standard k-e tended to overpre- 
dict secondary flow zones; the modified k-e model was postulated to provide 
better results. 

The standard k-e model tended to overpredict direct dynamic stiffness 
when compared with the modified k-e model; however, there was no change in 
quadrature dynamic stiffness. 

The M-shaped effective viscosity profiles suggest significant production 
of turbulence near the rotor at about 25 percent of the clearance and a lesser 
peak near the stator at about 75 percent of the clearance. The effective vis- 
cosity and peaks were maximum near the maximum clearance and diminished in the 
convergence zone, suggesting the possibility of relaminarization or superlami- 
nar flow. With viscous forces reduced in the minimum clearance region and 
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large in the turbulent flowfield outside this region, a physical model for sec- 
ondary flows could be developed. 

Without secondary flows the Richardson parameter was usually negative and 
decreased from the rotor to the stator; with secondary flows it first decreased 
and then increased, becoming positive near the stator and within the zone of 
secondary flow. A positive Richardson parameter implies that turbulent damping 
exceeds production, suggesting that designs where dissipation exceeds produc- 
tion near the stator are stable. As the rotor approaches the stator, an unsta- 
ble rotor can engender secondary flow patterns that quench the instability; 
this usually implies that the Richardson parameter near the stator changes from 
negative to positive and provides an explanation of a common class of turboma- 
chine oscillations cited by Childs. 
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(B) PASSAGE WITH ROTATING SHAFT. 

FIGURE 1. - SKETCHES OF TURBULENCE REDISTRIBUTION. 
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(B) QUADRATURE. 

FIGURE - DYNAMIC STIFFNESS FORCES FOR ±1:1 AND ±3:1 
INJECTION ANGLES WITH AND AGAINST ROTATION AS A FUNC- 
TION OF ROTOR SPEED, SYNCHRONOUS PERTURBATION, AND 
PRANDTL TURBULENCE MODEL. WITH DYNAMIC ECCENTRICITY/ 
CLEARANCE RATIO A/C OF 0.8. (WORKING FLUID, HALON 
(CBrF 3 ); PRANDTL MIXING LENGTH TURBULENCE MODEL ONLY.) 
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(A) ROTOR SPEED, 2344 RPM; INJECTION WITH ROTATION (1:1). 

FIGURE 5 - QUALITATIVE AXIAL AND CIRCUMFERENTIAL FLOW VELOCITIES AT FIVE AXIAL CROSS SECTIONS 
(L/A, L/2 - E. L/2, L/2 + E, 3L/4) WITH A/c = 0.8. (PRANDTL MIXING LENGTH TURBULENCE MODEL.) 

+ INDICATES w > 0. 
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(B) ROTOR SPEED, 234*1 RPM; INJECTION AGAINST ROTATION (-1:1). 


FIGURE 5. - CONTINUED. 
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(C) ROTOR SPEED, 5085 RPK; INJECTION WITH ROTATION (1:1). 
FI6URE 5. - CONTINUED. 
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(D) ROTOR SPEED, 5085 RPM; INJECTION AGAINST ROTATION (-1:1). 
FIGURE 5. - CONCLUDED. 
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FIGURE 6. - SEAL CONFIGURATION AND PARAMETERS FOR 
SWIRL FLOWS. 


KINEMATIC VISCOSITY, M z /s 


MODEL 


PRANDTL MIXING LENGTH 


MODIFIED k-£ 

STANDARD k-£ 

I.AxlO -1 * 




(A) 90° TO MAXIMUM CLEARANCE. (B) APPROACHING MAXIMUM CLEARANCE. 



RADIAL POSITION, r/c 
(C) APPROACHING MINIMUM CLEARANCE. 

FIGURE 7. - EFFECTIVE VISCOSITY AS A FUNCTION OF RADIAL POSITION AT THREE 
CIRCUMFERENTIAL POSITIONS. 
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NOI1F 1 ( NFAR ROTOR ) — i 



(A) PRANDTL MIXING-LENGTH MODEL. 



(B) STANDARD k-£ MODEL. 



(C) MODIFIED k-£ MODEL. 


FIGURE 8. - ZONES OF SECONDARY FLOWS AND CIRCUMFERENTIAL VELOCITIES NEAR ROTOR AND STATOR 
FOR THREE TURBULENCE MODELS. ROTOR SPEED. 5085 RPM. NO PRESWIRL. AND NO INJECTION. 
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L/4 SEAL EXIT 



NODE 6 (NEAR STATOR) 





(B) ROTOR SPEED, 3252 RPM; PRESWIRL AGAINST ROTATION. 

FIGURE 10. - QUALITATIVE CIRCUMFERENTIAL VELOCITIES AND VARIATIONS OF RICHARDSON PARAMETER 
WITHIN SEAL CROSS SECTION FOR SYNCHRONOUS PERTURBATION, A/c = 0.24, AND MODIFIED k-£ MODEL. 
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(D) ROTOR SPEED, 5085 RPM; PRESWIRL AGAINST ROTATION. 


FIGURE 10. - CONCLUDED. 
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ROTOR SPEED, RAD/s 


(A) DIRECT (K d ). 


(B) QUADRATURE (K Q ). 


FIGURE 11. - DYNAMIC STIFFNESS FORCES USING STANDARD AND MODIFIED k-£ TURBU- 
LENCE MODELS AS A FUNCTION OF ROTOR SPEED, SYNCHRONOUS PERTURBATION, AND PRE- 
SWIRL AGAINST ROTATION. £ = — . ( V PRESWIRL /R R U R = ±(1 - r/R $ ) 1/3 . ) 




(A) DIRECT (K d ). (B) QUADRATURE (K q ). 

FIGURE 12. - DYNAMIC STIFFNESS FORCES USING MODIFIED k-£ TURBULENCE MODEL AS 
A FUNCTION OF ROTOR SPEED, SYNCHRONOUS PERTURBATION, AND PRESWIRL AGAINST AND 
WITH ROTATION. (V pRESWIRL /R R U R = ±(1 - r/R s ) 1/3 .) 
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